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Abstract 

We study the drift of suspended micro-particles in a viscous liquid pumped back and forth through a 
periodic lattice of pores (drift ratchet). In order to explain the particle drift observed in such an experi- 
ment, we present an one-dimensional deterministic model of Stokes' drag. We show that the stability of 
oscillations of particle is related to their amplitude. Under appropriate conditions, particles may drift and 
two mechanisms of transport are pointed out. The first one is due to an spatio-temporal synchronization 
between the fluid and particle motions. As results the velocity is locked by the ratio of the space periodicity 
over the time periodicity. The direction of the transport may switch by tuning the parameters. Noteworthy, 
its emergence is related to a lattice of 2-periodic orbits but not necessary to chaotic dynamics. In this point 
of view it appears as a generalization of the Stokes' drift. The second mechanism is due to an intermittent 
bifurcation and leads to a slow transport composed by long time oscillations following by a relative short 
transport to the next pore. Both steps repeat in a quasi-periodic manner. The direction of this last transport 
is strongly dependent on the pore geometry. 
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I. INTRODUCTION 



sec:intro 

The transport of micro-particles through pores in a viscous fluid in absence of mean force gradient 
finds its motivation in many biological applications as the molecular motor or molecular pump. In 
the last decade, the literature shows that a periodical pore lattice without the symmetry x — > —x 
can lead to the so-called ratchet effect allowing an transport in one direction x or —x. A review 
can be found in IfTOll . In order to better understand the transport of particles in a pore lattice, the 
Max-Planck Institute in Halle (Germany) has build macroporous silicon wafer which is connected 
at both ends to basins [fT3l . The basins and the pores are filled with liquid and micrometer- sized 
particles pumped back and forth. Thus, a backward and forward motion of liquid occurs draging 
the particles. The experiment shows the existence of an effective transport in a certain range of 
parameter values. By tuning them, the direction of the effective transport may change. We focus 
on a simple deterministic model displaying these phenomena but not necessary for the experiment 
parameter values. 

The drift of suspended particles in a viscous fluid under an oscillating force of zero average is 
often explained by the well-known Stokes' drift [|29l . It is a deterministic mechanism involving 
the existence of a propagating wave of the fluid velocity field. The particle may drift by synchro- 
nization with the phase velocity of the traveling wave by a so-called phase locking phenomenon 
11231 . The drift velocity is the phase velocity. One may use the analogy with a surfer on a wave: 
the particle rides the traveling wave ifTTTl. This drift remains when noise is added, however the 
efficiency of the Stokes' drift decreases [[TTfl. In contrast for other systems the presence of noise 
plays a crucial role for the transport of inertial particles. It is noteworthy that the drift dynamic is 
during a time range close to the Stokes' drift: the drift velocity corresponds to the ratio between 
the spatial periodicity and the temporal periodicity of the forcing. Such a transport arises as well in 
Brownian particles [f3j[T6l as for electric current in a SQUID device [32]. Note, that this drift may 
exist without the presence of traveling wave of the force field. More surprising, if the average of 
the force is non-zero the drift may be in the opposite direction, it is the absolute negative mobility 

ma. 

For the transport of milli- and micro-particles in a viscous fluid, the role of the thermal noise is 
debatable. The fluid motion without particle is described by the Stokes equations (low Reynolds 
number) leading to an oscillating creeping motion. If the noise and the inertia are negligible the 
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particles advect with the fluid getting an oscillating motion [15J. Therefore, in order to explain 
particle transport, another phenomenon has to be take into account: either thermal noise or iner- 
tia, or both. Kettner et al. [[TBI neglect the inertia and consider thermal fluctuations. Stochastic 
simulations display transport motions and also the direction may change by varying the frequency 
of pressure oscillations. However, even the micrometric size of particles the inertia may be larger 
than the thermal noise. In this paper, we propose a complementary approach where the noise is 
neglected while the small inertia is taken into account because of a finite mass, then finite drag. 
Thus, it is a deterministic approach of inertial particle in a viscous fluid. The question of transport 
solutions is not trivial since the sinusoidal pumping leads to a flow without the presence of a trav- 
eling wave and then the system is different from the Stokes' drift. 

In order to clearly identify the mechanisms of the drift, we focus on the simple case where the 
competition between the flow drag and the particle inertia governs the motion. Thus, we propose 
an one-dimensional deterministic problem. The inertia in deterministic ratchet problem showed 
chaotic dynamics with possible transport and reversal current |3j. Mateos |[T8l shows that the 
current reversal for a biased oscillated force corresponds to a transition from a chaotic to a pe- 
riodic regime or chaotic crisis flU. It is a phase synchronization of chaotic oscillators described 
in the ratchet framework by lfT2l l20l ETl . This phenomenon is generic when a periodic external 
forcing is applied. It is observed in several experiments ranging. The resulting transport solution 
was thoroughly studied in ll28l in deterministic for unbiased force. They found transport solutions 
with locked velocity c = n/m with n and m integers. These solutions may exist for symmetric 
ratchets where transport are possible in both direction. A current reversal phenomenon is obtained 
by applying a non-zero biais because of different existence domains of the transport to the left or 
to the right. The emergence of the synchronized transport appears after evolving through a period- 
doubling route to chaos: attractor crisis [0. Another studies shows the same mechanism of the 
emergence of synchronized dynamics as OTI . However, it seems that chaotic transitions is not 
required for the existence of synchronized transport rather the coexistence of periodic orbits O. 
As for ll28ll . we will study the transport solutions for the symmetric case in order to point out the 
role of pore asymmetry underlying the ratchet effect. The synchronized transport solutions are 
tracked in the comoving frame using path-following method. Contrary to a direct time integration 
as it is carried out in the literature, one gets the stability of the dynamics the kind of transition and 
the domain of existence. For instance it is possible to observe unstable branches even if a stable 
branch attracts all the dynamics. This approach allows us to give a new interpretation of synchro- 
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nized transport solution as a generalized Stokes' drift related to 2-periodic orbits. Moreover, we 
are able to interpret the different transitions found in literature. Nevertheless the parameter domain 
for which the synchronized transport solution occurs is too small for micro-sized particles. For 
large drags there is still transport solution resulting from an intermittent bifurcation. This transport 
is slow and in contrast to the previous case is strongly related to the asymmetry. This transport 
does not seem appear in the literature. 

H. MODEL 

sec:mod 



A. Stokes approximation 

The Stokes equation is a model of viscous fluid for which the viscous dissipation is compen- 
sated with the pressure force work. Thus, it is a quasi-static fluid motion known as the creeping 
flow. The main parameter which is relevant for such an approximation is the Reynolds number 
which measures the mean fluid velocity scaled by the cinematic velocity. In order to determine the 
Reynolds number, let us give the magnitude order of the physical parameters of the experiment in 
Table [ij The fluid chosen in the water of density p = W 3 kg/m 3 and p = 10~ 3 Pa.s. The time 
is scaled by the period of the pumping and varies between [0.1; l]ms. The length is scaled by the 
pore length about L = 8.4/im. Finally, the velocity is scaled by the spatial mean velocity taken 
at the time where the flux pumping is maximal. One calls it the characteristic velocity U m . Using 
these parameters to scale the Navier-Stokes equation, one obtains: 

P U m Ldv pU 2 m T pTL UT 

ttt H (v.V)v = Vp+——Av (1) 

p at p p L 

(71) 

Re l — + Re 2 {v.V)v = -C x Vp + C 2 A 2 v (2) 
at 

The first two pre-factors Re\ and Re 2 are two possible expressions of the Reynolds number. Ac- 
cording to the Table [j] they are about 10~ 2 while both dimensionless parameters G\ and C 2 are 
about 1. Then, it is relevant to neglect the terms on the left side and one obtains the quasi- static 
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/i : Dynamic viscosity of the water W~ 3 Pa.s 

r p : Particle radius 1.8 • 10~ 6 m 

po : water density 10 3 kg/m 3 

p p : particle density ~ p 

U m : Characteristic velocity of the fluid (> 0) 2.1 • 1CT 2 - 1CT 1 m/s 

T : Period of the pumping 0.01 10 _4 s 

L : pore length 8.4 • 10 _6 m 

r : friction constant per unit of mass 10 6 s _1 

TABLE I: Dimensional physical parameters of the drift ratchet experiment. 

tab:expval 

approximation of the linear Stokes equation with no- slip boundary 

Av (r,t) = Vp(r,t), eq : stokes (3) 
v (r, t) = 0, on the pore boundary 

where Vq and po are, respectively, the velocity and the pressure of the fluid without particles. 
Because of the fluid is in equilibrium at each time t, the sinusoidal variation of the pressure implies 
a sinusoidal temporal dependence of the velocity field vo : 

Vo(r, t) = Uo(r) sin(27rt), eq : vxt (4) 

where w (r) is a velocity field obtained by applying a constant pressure difference back and forth 
of pore lattice. The lattice being periodic, then the function u (v) gets its periodicity. The adi- 
mensional characteristic velocity u m is the mean velocity u (r) in one pore. (All dimensionless 
parameters are noted using small letters). 

The simulation of the full non-linear Navier-Stokes equation with these parameters corroborates 
the relevance of the Stockes' approximation (Eq.|3]). 

B. Governing equation of the particle motion 

The velocity and acceleration of the particle motion being smaller than those of the flow motion, 
the Stokes approximation Eq.[3]remains valid. The flow motion is described by the Stokes equation 
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but only the boundary conditions of Eq. [3]have to be added in order to fit the no-slip condition on 
the particle surface: 

Avj = Vp/, eq : stokespart (5) 
Vf(f,t) = 0, on the pore boundary, eq: stokespart2 (6) 
Vf(r,t) = v(r,t), on the pore surface, eq : stokespart3 (7) 

where Vf(f,t) is the fluid velocity with the presence of the particles and v(r,t) is the particle 
velocity at the point r and the time t. At each time t the stress tensor on the particle surface 
may be estimated using boundary equations (see e.g. [@]|24]|). By adding the Newton's and the 
momentum equations of the solid particle, the particle motion can be completely described. The 
motion can be quite complicated if the particle does not stay on the symmetry x axis of the pore 
because, on one hand, the particle may rotate and, on another hand, the interaction between the 
particles and the wall has to be taken into account [5 |. Such a study is out of the scope of the paper 
since we aim at finding the most simple transport mechanisms in a viscous fluid. 
Therefore, we focus on the axisymmetric case: the particle is centered on the axis and does not 
rotate. With this assumption, the problem becomes ID. Let us note s(t) the position x of the 
particle center at the time t. Each point of the particle has the same velocity note v(t): 

m p d ^ = F(s(t),t), (8) 

where m p is the particle mass and F(s(t),t) the drag force on the particle. This last one is de- 
termined by the system [5j [6] and |7J Now one assumes that the drag force F is determined by the 
Stokes' drag of a spherical particle in an uniform flow, i.e. the drag force F is proportional to the 
relative velocity of the particle. Then 

F(s(t),t) = -0 (v(t) - v (s(t),t)) eq : linear - drag (9) 

The Eq. [9] is verified only by assuming: firstly, the pore boundaries are neglected in the drag 
force estimation and secondly, the unperturbated v flow velocity is almost uniform in the volume 
occupied by the particle. Therefore these last assumptions yield only for small particles about 0.1 
fxm in the framework of the experiment since the narrow region of the pore is about 2.5/im. 
In this case the particle motion is one-dimensional model described by the ODE 

^ = -^(v(t)-v (s(t),t)) (10) 
at m p 

= 7 (v(t) -v (s(t),t)) eq:ldim (11) 
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with 7 = The relevant parameter is the dimensionless drag 7 defined by 



7=^ = ^ (12) 
Introducing the mean particle density p p it comes: 

7 = |^ 03) 

According to the value of Table |I} we get the approximation of 7 ~ 10-1000. However, the lower 
value of 10 is obtained for particle radius about 5-6 jjm. Then, the Stokes' force can not anymore 
be determined by the linear relation Eq.[9]if the pore radius is still about 8 — 9/im. However, in the 
paper we do not restrict our study to the experiment framework, for instance the pore width may 
be larger in order that the Eq. [9] remains relevant. The main goal of this analysis is to answer the 
open questions, is it possible to transport passive particles in a slow viscous flow with zero mean 
value in time and space? And if the answer is positive, which mechanisms are responsible? In 
this way, one does not aim at computing the exact velocity field u (x) from the experiment set-up 
but the velocity field will be rather described by three parameters, for instance the characteristic 
velocity u m , the velocity contrast a and the spatial asymmetry d. These parameters vary around 
values relevant for Stokes' flow and small particle inertia. 

Finally, one obtains a second order differential equation of the particle motion 

s + 7s = 71*0 (s, t), eq : ldode (14) 

and according to Eq.|4mo(s, t) can be decomposed 

Vo(s,t) = Uq(s) sin(27rt) eq : ldode (15) 

where uq(s) is 1 -periodic function, i.e. the pore periodicity in dimensionless problem. 

C. Dimensionless parameters 

As explained above, the drag may vary from values about 10 to infinite. The dimensionless 
characteristic velocity u m defined by 

u m = v (x)dx (16) 
Jo 
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is about 1 to 25 depending on the pressure difference on both sides of the pore network. 
Then one introduces the scaled relative variations of the velocity w (x) by: 

u (x) = u m (l + aw (x)) (17) 

where w is a periodic function and smooth enough, i.e. C 1 of zero mean value: 

Wq(x + 1) = Wo(x) eq : wOper (18) 
w (x)dx = 0. eq: wOO (19) 



One assumes that \w (x) | and a vary in the range [0, 1] so that a is related to the extrema of u 

a = — . (20) 

Therefore, the parameter a is called the velocity contrast. 

The particle motion is governed by the second order non-linear ordinary differential equation: 

s + 7s = 7« m (1 + awo(s)) sin(27rf), eq : ode (21) 

Such a equation admits an unique solution C 2 for a given position and velocity (sj, i>j) at a time U. 
One can show that the particle acceleration s and its velocity s are bound (see Appendix I). 

III. BASIC RESULTS AND NUMERICAL METHOD 



There is no analytical results of the non-linear Eq. 21 in the general case. In this section, one 



gives the behavior of the particle in some limit cases. The generic case is numerically solved. 



A. Basics 



1. Uniform flow: a <C 1. 



If a = the flow without particles vo(x,t) is spatially uniform. The ODE Eq. 21 becomes 
second order linear equation with a periodic force and a damping force: 

s + js — '-yu m sin(27rt), eq : odea = (22) 

The dynamics tends exponentially to an 1 -periodic solution with the amplitude 
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The amplitude of the oscillations is maximal for large drag (advection). The characteristic time of 
relaxation is scaled by the inverse of the drag I/7. 

If a 7^ 0, the invariance in the x direction is broken and the existence of 1-periodic solution 
is no more guaranteed even if a is close to zero. Such a scenario occurs in the Section [V] for an 
asymmetric pore profile. However, if the geometry of the pore has the parity left-right symmetry, 
then there is generically solutions which have this symmetry. Because there is two points of 
reflection symmetry -the minimum and the maximum of the flow v (x)- there is two 1-periodic 
solutions which move around each symmetry point (see next Section [TV]). The existence of both 
solutions plays a relevant role for the understanding of the transport dynamics. 



2. Large drag 7. 

In the experience the drag is at least equal to 100. The limit case 7 = +00 corresponds to a 
pure advection of the particle and leads to the first order ODE 

Soo = u m (1 + cuuoOoo)) sin(27rt). (24) 

The solutions of this equations are 1-periodic in time. Then, a drift is not possible if no constant 
force is present at the gravity. In the literrature it will be found as overdamped ratchet. In the 



Section V D we will study the nearly case of overdamped ratchet: 7 is finite but large. 

For large values of 7 but finite, the continuum of 1-periodic solutions disappears but it may 
still remain both 1-periodic solutions mentioned in the paragraph above. The existence of periodic 
solutions at large drag 7 will be discussed in the last section [V] and it depends strongly on the 
asymmetry of the flow. 

Then, for small values of parameters a, u m and large values of 7, one expects only 1-periodic 
motions when the flow profile is symmetric or slightly asymmetric. 

3. linear variation of vq(x) 

Here, we do not consider a pore lattice rather a pore with varying diameter so that the velocity 
v varies linearly 

x 

u (x) = u m (l + — ), eq : wOlin (25) 
a 



Prenrint - March 4. 201 3 



with a > and u m > 0. For a dynamic close to x = 0, for instance far from the singular point 



x = —d, the time integration of the ODE (21 ) shows a drift to the negative direction, i.e. to the 
decreasing values of v (x). The consequence for the particle motion in a pore is that for small 
enough oscillations, the particles tend to the region where v (x) is minimal, i.e. where the pore 
has the larger diameter. Then, the periodic solution close to the velocity minimum should be 
stable and the other one close to the maximum should be unstable. Of course, when the oscillation 
amplitudes become comparable to the pore length, this argument does not yield any more. 

B. Numerical method 

The above analysis shows that it may exist two 1 -periodic solutions if the profile is slightly 
asymmetric and the oscillation amplitude is small. Moreover, one does be stable. Therefore, 
this scenario does not allow the existence of a drift. We need to explore larger parameter values. 



The time integration of the Eq. 21 allows to find possible transport solution but each simulation 
depends on the initial values and there is no information on stability. Moreover, it is difficult to 
understand the mechanism of the transport. The continuation technics is a powerful tool to track 
periodic solutions in the parameters space. Moreover, one gets information about their stability. 
This approach allows to better understand the dynamics and to find transport solution. In some 
cases, the onset of existence of periodic solutions may explain a transport phenomenon. In order to 
perform the continuation task, one uses the freeware package AUTO [7J. It allows the continuation 
of periodic orbit and the detection of doubling period and fold bifurcation of a ODE's system. The 



ODE system (21 ) is transformed in a three dimensional system with (s, v(= s),t) as variables 



s = v 

v = 7 (u m [l + awo(s)] sin(27rt) — v ) (26) 
i = 1 

This last system can be formal written 

(s, v, t) = F(s, v , t) eq : ds (27) 

Each point of the bifurcations diagrams represents the norm |.| of a T-periodic state s. Its norm 
is defined as the L 2 -norm of the particle velocity: 

,t \ 1/2 

10 



| s || = I ( / v \t)dt ) . (28) 
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Then, the norm does not depend on the particle position and for instance two identical particle 
motions but shifted by an integer of pore lengths have the same norm. 

The path-following method requires to know the starting solution. The periodic motions for an 
uniform flow, i.e. a = 0, being known, one performs the continuation of the periodic solution 
branch expanding first by the parameter a. 



IV. PARITY-SYMMETRY OF PORES 



sec:sym 

If the geometry of pores lattice has the spatial parity symmetry x — > —x, an effective transport 
is not possible. However, this fact does not forbid the existence of transport solutions for a given 
initial condition. Then, the symmetry of the flow velocity implies the existence of a symmetric 
transport in the opposite direction. The motivation to consider a parity symmetry geometry lies 
in a better understanding of the different mechanisms resulting in a ratchet effect: dissymmetric 
flow and synchronization phenomenon. In this section, we show the existence of transport so- 
lutions related to a spontaneous symmetry breaking of symmetric periodic solutions. In the next 



Section V B the existence of transport solution for an asymmetric geometry is interpreted as a per- 
turbation of the symmetric case using continuation method. The asymmetry of the pore geometry 
ensures only that the mean transport measured over all initial conditions is generically non-zero. 



A. Equivariance of the ODE Eq. 27 



The parity symmetry means the parity of the incident velocity field at a given time t, i.e. 



V (-X) = Vq(x) <t==^ W (-X) = Wq(x) 



(29) 



Then, the function F of the system Eq. 27 is equivariant by the spatio-temporal transformation 



Sq : (s, v, t) — y (— s, — v , t + 1/2) eq : equiO 



(30) 



what means 



F(S {s,v,t)) = S {F(s,v,t)). 



(31) 



The consequence is that for each solution (s(t),v(t)) of Eq. 27 then the symmetric dynamic: 
(s-(t), = (— s(t + 1/2), — v(t + 1/2)) is solution too. Let us remark that because of the 
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time and spatial 1 -periodicities, the system is also equivariant by the transformation 



S m : (s — 1/2, v, t) — > (— (s — 1/2), — v,t + 1/2), eq : equidemi (32) 

which corresponds to a parity symmetry around the middle of the pore: x = 1/2. 

We say that a periodic solution is symmetric if it is invariant by one of these transformations S or 

S m . We denote by s the solution invariant by S . If < s > t is the mean position over the temporal 

period: 

1 <- T 



< s > t = - / s(t)dfc. (33) 







then the equivariance show that < s >t— 0. In this sense it is centered at zero. In the same 
way the symmetric solution s m invariant by S m is centered at the pore middle: < s m > t = 1/2. 
Obviously, the spatial periodicity implies that at each point x — k/2, k E Z there is a symmetric 
solution. But there are spatially translated solutions from s or s m . Then we focus only on both 
solutions s and s m . 

Because of the spatio-temporal symmetries, the solutions s and s m generically exist according to 



the concluding remarks of Section III A In order to find transport solutions, this parity symmetry 
of the dynamics needs to be broken. Then by varying the three parameters: the drag 7, the charac- 
teristic velocity u m and the velocity contrast a, we explore the bifurcation diagrams of the periodic 
solution branches. 

A sinusoidal dependence is chosen for the velocity field v (x) of the driving flow: 

v (x) = u m (1 + a cos(27rx)) , eq : uOsym (34) 

which corresponds for a = 0.65 roughly to the velocity profile on the axis in the micro-pump 
experiment in lfT3l . Thus, this value will often used for the different bifurcation diagrams. 



B. 1-periodic solutions 

The parity- symmetric periodic solutions s and s m get the periodicity of the forced oscillations, 
i.e. they are one-periodic. The 1-periodic symmetric solution is given for a = by the linear 



differential equation Eq. 22 and it is the starting solution used for the path-following. 

The first bifurcation diagram is performed with a large drag 7 = 100. Then, one expects an 

almost advective motion of the particle of small amplitudes. When a is non-zero the continuum 

12 

Prenrint - March 4. 201 3 



0.6 




-0.2 ,-' 



-0.3 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 

Time t 



FIG. 1: Time evolution of the 1-periodic symmetric solutions (red) sq and (black) s m for u m = 1, 7 = 100 



and a = 0.65 of the velocity field Vq(x) given by Eq. 34 Solid [Dashed] line indicates stable [unstable] 
solution. 

fig:symulsolt 

of periodic solutions breaks into the two solutions sq and s m as shows Fig. [T] The solutions so 
and s m are respectively oscillations around the maximum and the minimum of the velocity field. 



As expected in Section III A|the solution s is unstable and s m is stable. The continuation of these 



solutions shows that for any velocity contrast a there is no bifurcation and the stable solution s m 
remains stable. By decreasing the drag, there is no bifurcation too when a = 0.65. Therefore, all 
the dynamics are attracted by the only stable periodic solution s m . 

Now one varies the characteristic flow velocity u m and the velocity contrast and the drag are 
fixed: a = 0.65, 7 = 100. Both solutions exchange their stability (Fig. [2]). The solution s 
stabilizes at = 5.27. The Fig.[2]-b shows its time evolution during one period: the amplitude 
ends just beyond the position x = ±1/2, i.e., where the flow velocity is minimal. In contrast, 
the periodic solution s m destabilizes at u™ 1 = 4.84 for which the amplitude of oscillations ends 
close to x — 1, i.e. the maximum of the velocity flow (Fig. [2]-b). Moreover, by increasing u m , 
the solution stabilizes again at w™ 2 = 9.24 when the amplitude of oscillations are just beyond the 
minimum of velocity x — 1.5 or x — —0.5 (Fig.|2]-b). An analogous destabilization occurs for s 
at = 8.79: the boundary of the oscillations of the periodic solution s are close to ±1 (maxima 
of flow velocity). Therefore the stability of the one -periodic solutions seems to follow the rule: if 
the boundary of oscillations are close to the minimum of velocity vq (i.e. at x = l/2mo<i[l]) the 
solution is stable and if the boundaries of oscillations are close the maximum of the velocity v , 
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FIG. 2: Continuation of branches of 1 -periodic orbits (red line) sq and (black line) s m for the symmetric 



velocity profile Vq(x) (Eq. 34 1 by varying u m at 7 = 100 and a = 0.65. Red color corresponds to sq and 
black to s m . Symbol 'x' indicates a pitchfork bifurcation, (b) Time evolution of periodic branches at the 
pitchfork bifurcations of the panel (a): (plain red) so at u™ 1 = 5.27, (dashed red) sq at -u™ 2 = 9.24, (dashed 
black) s m at vfy =4.84 and at (plain black) u% =8.79 

fig:symgia65cu 

then the solution is unstable. 

The exchange of stability implies the emergence of a new branch which breaks the Sq reflection 
symmetry, then it is called the asymmetric branch s a . It is a pitchfork bifurcation of periodic 
orbit and then there are two solutions s a+ and s a _ which emerge from the critical point such 
s a + = So(s a -). Because they are symmetric, they have the same norm and then they do not 
appear twice in the bifurcation diagram Fig.[2]-a. Thus both solution branches are called s a . When 
the solution s m loses its stability at w™ 1 the branch s a emerges supercritically and connects the 
branch s at ify which gains stability. A similar scenario appears between the critical values 
and w™ 2 (Fig. [2]). The branches s a being stable, there is always a stable periodic solution for all 



values of u m . The time integration of the ODE Eq. 21 shows that the dynamics tend to the stable 
solution which forbids any transport solution. 



To conclude, the symmetric solutions s and s m may change their stability which depends 
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FIG. 3: Continuation of 1 -periodic orbits by varying the drag 7 from 100 to zero. The starting states at 
a = 0.65 are the circles in the bifurcation diagram Fig. [2^. Color code of the branches is the same as 
Fig. [2^. Symbol 'x' indicates a pitchfork bifurcation while a dot indicates a period doubling bifurcation. 

fig:symu5a65cg 

strongly on the amplitude of their oscillations. Roughly, the stability changes when the oscillation 
amplitude increases of a half pore length. At the onset of stability, a new branch s a of 1 -periodic 
solution bifurcates. This branch is no more invariant by the reflection symmetry S nor by S m . 
For a large drag (7 = 100) this asymmetric branch s a is stable and then there is always a stable 
periodic solution which attracts the dynamics. In the following, one studies the dynamics into 
smaller friction term 7. 




C. Doubling period cascade 

In this section, we show that the s a branch may be unstable and then a cascade of period 
doubling bifurcations occurs as it is generically exhibited by nonlinear oscillators with dissipation 
and additional time dependent force lf30ll . For that let us carry out the continuation of the three 
1-periodic solutions at a = 0.65 and u m = 5 (circles in Fig. [2^) by decreasing the drag. The Fig. [3] 
shows that the branch s a connects the branch s which gains stability at 7 m = 5.67. In contrast to 



other scenarios of Section IV B the branch s a becomes unstable in the range pyjx = 6.67; 71^ 
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FIG. 4: Bifurcated period-doubling branches emerging from (green line) s a in Fig. [3] Shown is (a) the 
bifurcation diagram and (b) the magnification near the period doubling bifurcations. For the one-periodic 
orbits the color code is the one used in Fig. [3] The 2-, 4- and 8-periodic orbits correspond respectively to 
the thin blue, cyan and sky-blue lines. Solution profiles at 7 = 10 indicated by circles are shown in Fig. [5] 

fig:symu5pd 

16.0] via two period doubling bifurcations. Both period-doubling bifurcations are connected by 
a branch of 2-periodic orbits (see panels (b) of Fig. [4]). This last kind of solution needs two 
temporal periods to come back to the same state (position and velocity). In particular, there is 
two different maxima (Fig. [5]). Due to the finite drag, there is a small delay for the direction 
change of the particle since the extrema are slightly after the times t = nrr/2 with n entire. 
The 2-periodic branch is stable close to the bifurcations points but loses its stability via a period 
doubling bifurcation at a slightly larger or smaller 7. Again a 4-periodic branch relies these points 
(74X, = 6.98 and 74^ = 13.8) and it loses its stability via period doubling at 7sl = 7.06 and gains 
its stability at 'jgR = 13.5. Then, there is a range for which the 4-periodic solution is unstable. The 
8-periodic branch is displayed in Fig.|4]-a and its period-doubling bifurcations are very close from 
the branches emergence (Fig. [4})). This cascade of period doubling continues without end and 
leads to chaotic dynamics. Moreover, the period doubling bifurcation of the branch 2 n+1 -periodic 
is always closer to the previous bifurcation point of the 2 n branch implying that there is a range 
of non-zero length where every periodic branch is unstable. Thus, there is a range for which all 
periodic orbits are unstable and letting appear quasi-periodic or chaotic motions. The poincare 
section in the plane (x,v) in Fig. [6] shows a strange attractor for 7 = 11. The visited domain is 
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FIG. 5: Time evolution of the l-(plain green line), 2- (blue dots), 4-(magenta 'x' symbols) periodic solutions 
for 7 = 10 at locations indicated by circles in Fig. [4] 

fig:snapd 

included in one pore length thus it is bound. The dynamics may be a quasi-periodic since there is 
often the case for systems with forced oscillations. This route to the chaos is widely described in 
deterministic inertia ratchet fl^[19l|27l[28]|. Note however that the route to chaos studied in these 
papers concerns already drifting solutions which are periodic in the comoving frame. 
Eventually, one found also 3— and 5— periodic solutions. These solutions may be stable but they 
lose their stability by the same period-doubling phenomenon. We will study such a dynamic for 
non symmetric geometry in Section fV} 

The complexity of the dynamics exemplifies the classical results relative to forced systems 
present in a plethora of nonlinear oscillations topics. Here we do not aim at studying chaotic 
dynamics since they do not necessary lead to a transport solution. 

D. Synchronized transport solutions 

In order to introduce the transport solutions, one tracks the 1 -periodic solutions by expanding 
the mean velocity u m at 7 = 10 and a = 0.65 (Fig. [7]). The branches sq and s m exchange their 
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FIG. 6: Strange attractor in the (cc, v) Poincare map of the dynamics at 7 = 11 for the parameters of the 
bifurcation diagram Fig. [4] obtained for 100 different initial conditions. Parameters: a = 0.65, u m = 5. 

fig:symu5quasiper 

stability via the asymmetric branch s a as in the bifurcation diagram Fig. [2} However this last 
branch s a loses its stability by period-doubling bifurcations for u m about 5 or 10. One focuses on 
the period doubling cascade when u m is about 10. The s a branch loses its stability at u 2 ^ = 8.17 
where a 2-periodic branch appears (Fig. [8]). This last branch has 2 saddle-nodes (at u^ 1 = 10.64 
and at u^ 2 = 8.534) and connects the branch s a at w^f = 11.18 which becomes stable. The birth 
at u 4 ^ = 8.217 of the 4-periodic branch is very narrow to the 2-periodic one and connects the 
2-periodic branch at uf^ sn closer than 5 • 10 _7 % to u^ 1 the upper turning point. Thus in the short 
range [w^f sn ; (not visible in the Fig. [8]) the 2-periodic branch is stable. A second 4-periodic 
branch emerging from the 2-periodic branch at w^f = 11.10 connects the 2-periodic branch at 
u^ sn close to its second saddle-node u™ 2 : uf^ sn is about 2 • 10 _8 % superior to u™ 2 . Thus, again 
there is a narrow range [m™ 2 ; uf^ sn ] of stability of the 2-periodic branch. It is noteworthy that 
the turning points of the 2 and 4 periodic branches end roughly at the same value. Indeed this 
phenomenon occurs also for the 8-periodic solutions (not shown). These branches emerge very 
close to the 4-periodic branches and therefore, expect for islands of stability close to saddle-nodes, 
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FIG. 7: Continuation of 1-periodic orbits branches spanned by the characteristic velocity u m for a = 0.65 



and 7 = 10. The velocity profile is given by (Eq. 34). The same color code and the same bifurcation 
symbols as in Fig.[3]are used. 

fig:symcu 

a window of instability of periodic solutions occurs for approximatively 8.3 < u m < 10.7. The 
upper limit of this range is due to the branch s m which is stable after a turning point (see Fig. [7]). 
In this region we are able to observe quasi-periodic motion as shown in the previous section. 
However, the simulation shows a regular and stable transport solution at u m = 9 too (Fig. [9]). After 
one time period the particle moves of one spatial period. Thus the transport solution s t (t) seems 
to be the sum of an uniform translation at the velocity 1 plus an 1-periodic solution s p (t): 

St(t) = s p (t) +t eq: twsum (35) 

with 

Vt, s p (t + l) = s p (t) eq : sper (36) 

In other terms s t verifies 

VA; G Z, Vt G R : s t (t + k) — s t (t) + k eq: twper (37) 
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FIG. 8: Magnification in the range u m S [8; 11] of the bifurcation diagram Fig. [^(parameters: a = 0.65, 
7 = 10). Shown are the (dark blue) 2- and (magenta) 4-periodic branches emerging from the (green lines) 
s a branch. Unstable branch is represented by dashed line otherwise the line is plain. The color code and the 
bifurcation symbols are the same as used in Fig.|4| 

fig:symglOcu 

This transport solution will be called synchronized transport solution and 1 -periodic transport 
solution too. These phenomenon is known as phase synchronization or phase locking effect DUEl 
[18l|26l|28l. In point of view of transport of particle, it may be interpreted as a modulated Stokes's 
drift which occurs for a traveling wave potential. In our case, only the periodicity of the potential is 
required and so it generalizes this transport kind. In order to track this solution and its properties, 
for instance its existence domain and its stability, one carries out the path-following of s p , i.e. the 



transport solution s t in the comoving frame. If s t is solution of the ODE (21 ) then s p is solution of 
the new ODE: 

s P + 7S P = 7 [u m (l + aw (s p (t) + £)) sin(27rt) - 1] eq : tw (38) 



The presence of the constant biais —1 on the right side of the Eq. 38 recalls the tilted force in 
ratchet problem [25 J. This force tends to move the particle to the right what corresponds to periodic 
solution s t in the laboratory referential. This problem is analogous of current reversal phenomenon 
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FIG. 9: Time evolution of the synchronized transport solution st at u m = 9. The remaining parameters are 
as in Fig.[|(7 = 10, a = 0.65). 

fig:symtw 

in tilted inertial ratchet (see e.g. [19]). Here we seek periodic solutions in the comoving frame so 
one aims at stopping the "natural" transport. 

The continuation of 1 -periodic solution s p solution of Eq. [38] expanded by u m shows a existence 
domain in a eight-shape delimited by the range w^ nl = 8.65 and w^ n2 = 10.62 corresponding to 
two saddle-nodes (Fig. [TO]). Then there is always two different solutions and at most one solution 
is stable. On the left, the first range of stability starts from the turning point and ends at the period 
doubling bifurcation at u 2 ^ = 8.676. The second range of stability is delimited by two period 
doubling bifurcations at u 2 ^ 11 = 8.97 and u™ 2 = 9.37. Finally, the third range is very narrow 



(not visible in Fig. 10) between u 2 ^ = 10.61 and the rightmost saddle-node a™ 2 . In these stability 
regions the transport solution is stable. Note, that the existence of a transport solution s t+ in the 
direction to x positive implies the existence a transport solution s t _ in the negative direction given 
by: 

s t _(t) = S (s t+ (t)) = -s t+ (t + 1/2) (39) 
st-(t) = s p (t + 1/2) - 1/2 - t eq: tw - sum (40) 

If for an initial condition (xi,Vi,U) the dynamics tends to the s t+ solution then for the initial 
condition (— x%, —Vi,U + 1/2) the dynamics tends to the s t _. Therefore, there is four 1-periodic 
transport solutions as mentioned in [|28l . In the three stability ranges, there is two stable transport 
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FIG. 10: Continuation of the periodic orbits s p (Eq. 38 1 associated to the synchronized transport solutions 



St with drift velocity c = ±1. Black dots are period doubling bifurcations. The circles indicate the locations 



of solutions shown in Fig. 1 1 Blue lines are 2-periodic transport solutions. The remaining parameters are 



as in Fig. §(7 = 10, a = 0.65). 

fig:symctw 

solutions in opposite directions that attract all the dynamics. 
Let us define the transport velocity c by 

x(t)-x(t ) 

c = lim eq : c (41) 

t^oo t — to 

if this limit exists. Otherwise, we introduce the transport velocity at the time t: 

x(t) — x(t ) . . 

c(t) = y, t _^ ° J eq : c(t) (42) 

According to the Eqs. ([35]) and ([40]) it is clear that the transport velocity c is locked to c = ±1 for 



all solutions of the branch displayed in Fig. 10 



In the instability domains a period doubling cascade occurs for the solution s p similar to the 
figures [4]-a or [8] Fig.[10]displays the 2-periodic branches which connect the four doubling period 
bifurcations. An example of the resulting transport in laboratory referential corresponding to the 



2-periodic solution s p is shown in Fig. 1 1 Two periods are required in order to shift exactly the 
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FIG. 1 1 : Time evolution of the synchronized transport solution st at location indicated by circles in Fig. 10 
Plain black line corresponds to the solution at u m = 9 and the blue + correspond to the 2-periodic transport 
solution st at u m = 8.93. The remaining parameters are as in Fig.[8](7 = 10, a = 0.65). 

fig:symtwpd 

particle of two pores. Thus, the transport velocity remains equal to one. In a more general way, 
all T-periodic solutions of the Fig. [10] (in co-moving frame) implies the existence of a transport 
solution with |c| = 1. The period doubling cascade implies the existence of a parameters range 
where all periodic orbits s p are unstable and then the transport solution s t too. In this region, 
the dynamics switch alternatively between transport solution to the right direction and transport 
to the left. An example in given in Fig. 12 for u m = 8.8 with two different initial conditions. 
During a few tens of periods the dynamics are close to the transport solution with c = 1 or — 1 
and the change of the direction needs only one or two time periods. We do not observe dynamics 
close to periodic or quasi-periodic motion (laboratory frame) which are unstable too. The discrete 
dynamics at every entire time n is similar to a random walk. The velocity transport remains 
close to zero as shown the simulation for u m = 8.8 with different initial conditions (Fig. 12). 
Therefore there is a brutal transient of the dynamics from a transport at c = ±1 (stable periodic s p 
solution) to a zero mean drift c = (unstable bound solutions s p ). Thus the discontinuity of the 
drift velocity is related with a period doubling and it is known as a chaotic transition or a crisis 

mm. 

The existence domain of synchronized transport solution is represented in the plane (u m , 7) in 
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FIG. 12: (a) Particle dynamics in the presence of two unstable opposite transport solutions at u m = 8.8 for 
two different initial conditions. The remaining parameters are the same as in Fig.[8](7 = 10, a = 0.65). (b) 
Convergence of the transport velocity of the dynamics of panel (a). 

fig.symtwt 



Fig. 13 The stable region of transport splits into many branches for 7 than 8 (not shown) and 
it is very reminiscent of the figures in [|28l . Indeed in this latter paper, the figures result from 
a time integration so that only transport solution having attractive manifold appear. It is the 
advantage of the continuation approach which allows displaying all branches even unstable. The 
transport solution does not exist for 7 superior to 13 which is still a small value for the experiment 
described in the introduction. Beyond this existence domain there is transport events resulting in 
a complex dynamics. Because of the parity symmetry the mean transport is zero so that we study 
this dynamics for an asymmetric geometry. 



The bifurcation diagram Fig. 10 corroborates that synchronized transport is "born out of 
tangent bifurcation" 11281 . Usually, the end of periodic transport is interpreted in literature as a 
way of crisis flUEEHEHl. ^ ne contmua ti° n of the transport branch shows even this crisis that the 
solution still exists but it is unstable. Then, the synchronized transport emerges and ends at a 
tangent bifurcation. This bifurcation being generic, it is surprising that such a periodic transport 
emerges from bound chaotic dynamics such shown in Fig. [6] In fact, the 2-periodic branch plays 



a key role in the emergence of the synchronized transport solution. Let us represent in Fig. 14 



the 2-periodic solution at the emergence of the transport solution u m = (the right turning 
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FIG. 13: Existence domain of the synchronized transport velocity with c = ±1 demarcated by the bold line. 
Both subdomains (fine lines) bound the unstable solution and the remaining domain between bold and fine 
lines corresponds to stable solution. 

fig:twrange 

point of the Fig. [9]). Among the three 2-periodic orbits, one takes the orbit presenting the largest 
shift between both maxima. This s 2 solution is closed from the turning point and belongs the 
branch between both turning points (a™ 1 , «™ 2 ) of Fig. [9] In the same graph, we represent the 
s^p solution obtained from s*p by a space and a time period shifts, i.e. s 2 (t) = so(t + 1) + 1. 



This solution exists due to the discrete spatial and temporal translations symmetries. Fig. 14 
points out that trajectory of two consecutive orbits (t) and s 2 (t) are close around the time 
t = 0.25 + 2k, k e Z. At this time the velocity (the slope) is maximal and the particle position is 
roughly in the middle. Because the slope of both curves are almost equal, both trajectories in the 
3-dimensional phase space (x, v, t) are close together. If the 2-periodic solutions is not stable, the 
dynamics in the vicinity of s 2 °^ may switch to the s 2 at t about 0.25. The representation of the 
transport solution corroborates this fact. The synchronized transport solution s t is close to for 
t between [0; 0.25]. At t about 0.25, the values of s t are close to both curves s 2 an d s 2 . Then, s t 
moves slowly away from and brings closer to s 2 . The trajectories of periodic orbits diverge 
after passing the narrow region of pore at x — 0. The s t solution follows the s 2 orbit and stays 

close to s 2 till t ~ 1 + 0.25. Because of the spatio-temporal symmetry a similar scenario allows 

(2) 

the trajectory to switch to s 2 . The "bridge" between two consecutive 2-periodic solutions being 
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FIG. 14: Time evolutions of (dots) synchronized transport solution and (plain and dashed lines) and 2- 
periodic solutions S2 at u m = in Fig. [8] Each solution s' 2 (dashed line) is deduced from the solution S2 
by spatial and temporal translation: s' 2 (t) = S2(t + 1) + 1. Parameters 7 = 10, a = 0.65. 

fig:gluing 

spaced of length one, thus one obtains a transport velocity of exactly one unit. 

Such a key role of the array of 2-periodic orbits in the transport solution existence is observed 



for u 



tsnl 



< U m < U 



tsn2 



The transport solution is close to the 2-periodic solution which is along 



the branch joining both turning points u s £ and uf% in bifurcation diagram Fig. |8j Moreover 
both turning points of this last branch corresponds roughly to the domain existence boundaries of 
synchronizes transport solution. This fact may be explained by plotting the solution of the unique 
2-period branch when u m < u m n2 or u m > u m nl . The difference between maxima or minima 



of this 2-periodic orbit is small (Fig. 15) contrary to the solution displayed in Fig. 14 Indeed 
this difference has to approach one in order to obtain a neighborhood between consecutive orbits. 
Therefore, the existence of a common neighborhood of two consecutive 2-periodic orbits seems 
to be the key of the synchronized transport existence. Thus, the whole bound chaotic attractor is 
not relevant. This remark is in agreement with Barbi and Salerno [2] who point out phase locking 
phenomenon without chaotic transition. 
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FIG. 15: Time evolutions of 2-periodic orbits at u m = 10.7 (circle in Fig. [8)7 = 10, a = 0.65). The orbit 
is shifted from sq such s\(t) = so(t + 1) + 1. 

fig:noglu 

In a more general way, any synchronized transport solution has the discrete velocity drift 
c q = l/q,q G Z may be emerged from of an array of m-periodic orbits with m > q + 1. Note 
that for larger characteristic value of u m the velocity may be rational c pq — p/q with p, q E Z. 
However, a characteristic velocity much larger than 10 is not in agreement with the Stokes 
approximation so that one considers only velocity c q inferior to one. Such a stable dynamics 
for c = 1/2 was found for u m = 9 and 7 = 13.153. The continuation of this transport solution 



by varying the velocity u m is shown in Fig. 16 1. The temporal evolution of the stable and an 



unstable transport solution corroborates the relation s t (x,t + 2p) = s t (x + p,t) in Fig. 17 In 
the same figure, we represent 3-periodic solutions which are translated of one spatial period and 
two temporal periods, i.e. the solutions s(x + p, t + 2p), p E Z. The unstable transport solution 
is very close to the periodic solution during two periods and then switch to another orbits. As for 
c = 1, two consecutive orbits are close where the velocity is maximal (middle of the pore) and the 
acceleration is small (linear variation). It is in this latter region that the switching occurs. Note 
that the stable solution is not very close to the orbit, nevertheless the switching occurs in the same 
region pointing out the same mechanism. Moreover, the existence domain of transport solution 
(Fig. [16^) and the one of 3-periodic solutions which have a large maxima difference (Fig.[T6p) are 
related. The existence domain is roughly the range between two saddle nodes of the 3-periodic 
orbit. Beyond the right saddle-node, both 3-periodic branches display orbits which have a small 
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(a) (b) 

FIG. 16: (a) Branch of transport solutions with c = ±1/2 with u m as parameter bifurcation, the other 
parameters being fixed: 7 = 13.153, a = 0, 65. The velocity field uq symmetric, (b) Branch of 3-periodic 
orbit with the same bifurcation parameter u m and the same fixed parameters. 
Dashed lines indicate unstable solutions in both bifurcation diagrams. 

ftg.cl/2 



difference between maxima resulting in a fairly spaced orbits array as in Fig. [T3J That explains 
that there is no transport solution. 

Finally, note that the stability domain of the transport solution is quasi-discrete so that one 
observes this solution for a specific parameter set like a resonance phenomenon. In the domain 
where the transport solution is unstable an intermittent dynamics may occur. In this latter case, a 
direct time integration of the governing does not display clearly a transport solution. 



In conclusion, the symmetric so and s m 1 -periodic solutions exchange their stability via a pitch- 
fork bifurcation. For drag about 10 the asymmetric solution s a presents a cascade of period dou- 
bling leading to chaotic but almost bound dynamics. In this parameter domain a synchronized 
transport known as a phase locked transport is found. The velocity of the transport is locked to 
plus or minus one. This synchronized transport may be stable and becomes unstable via periodic 
doubling. The transport dynamics then switch from c = +1 to c = — 1. Tangent bifurcations 
bound its domain of existence. This tangent bifurcation occurs in a chaotic region however the 
emergence of transport solution is favored simply by the array of 2-periodic orbits obtaining by 
the spatial and temporal symmetry translations. Two consecutive orbits having locally a common 
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FIG. 17: Temporal evolution of solutions represented in Fig. 16 by a circle at u m = 9. The plain and 
dashed black lines are 3-periodic solutions translated in space and in time. The trajectory with dashed line 
szd is deduced from ss p (plain black line) by: ssd{t) = + 2) + 1. The red plain line and the dotted 
trajectories are the stable and an unstable transport solution, respectively. 

fig:cl/2t 

vicinity. This vicinity happens when there is a difference between extrema of the periodic orbit 
is about one. Therefore such transport solution is characteristic of period doubling phenomenon 
in the presence of a spatial and temporal periodicities. The asymmetry of the transport solution is 
due to a spontaneous symmetry breaking of the 2-periodic solutions. Another transport solution 
with c = ±1/2 was found too. Its existence is related to the 3-periodic orbits in the same way as 
transport solution with c = ±1 is related to 2-periodic orbits. 

Beyond and close to the existence onset of synchronized transport an intermittent transport is ob- 
served. These latter dynamics are studied in the asymmetric velocity flow for which a preferential 
direction over all initial conditions is possible. 

V. ASYMMETRIC PORE 



sec:asym 
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FIG. 18: Velocity profiles of the driving flow uq{x) for u m = 1 and a = +0.65 (left panel a) and a 
(right panel b). The asymmetry d parameter takes the value —0.2, 0, 0.2 and 0.4. 

fig:velo-asym 

A. Introduction 



(b) 



-0.65 



In this section, one breaks the parity- symmetry of the flow, i.e., wo(— %) 7^ uq(x) an d then the 
existence of a transport solution to one direction does not imply the existence of a symmetrical 
transport solution in an opposite direction. Thus, an effective transport of particles is possible: 
the mean of velocity transport on all initial conditions is non zero. In a first part, we discuss the 
survey of the synchronized transport. In the second part, we show the existence of intermittent 
transport. We distinguish a chaotic intermittent transport for 7 about 10 and quasi-periodic inter- 
mittent transport for large 7. 

In order to have a dissymmetry parameter d, an asymmetrical flow u$ is constructed by the con- 
catenation of two functions cos with different scales: 

1 



U [X 



U r 



1 + a 



COS 7T 



X 



1/2 + d 



IAx) + COS 7T 



X 



1/2 -d 



'l-h{x)) 



eq : voasym 
(43) 
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with 

x = x mod (1) (44) 

f 1 if x < 1/2 + d 
k{x) = { (45) 

I otherwise 

The parameter d of asymmetry varies in the range ] — 1/2, 1/2 [. The resulting velocity field is 



represented on the Fig. 18 For d = the velocity profile is symmetric as in Eq. 34 When d 
is positive the minimum of velocity field is shifted to the right: in the region x E [0; 1/2 + d) 
the velocity gradient is smaller and in the remaining region [1/2 + d; 1] the velocity gradient is 
larger. The case d negative is deduced from the case d positive by the transformation x — > —x. 
Thus, in particular, a drift solution for d > implies a drift solution in opposite direction for 



d < 0. Moreover, contrary to the symmetric case (Fig. 18), negative velocity contrast a leads to 
a new velocity profile. However, one can show that if s + (t) is solution for a velocity contrast a 
then S-(i) = —s(t + 1/2) + 1/2 + d is solution too for the opposition velocity contrast —a. For 
instance, if s + is a transport solution for a, then s_ is a transport direction for —a in opposite 
direction. Therefore, one will study the influence of the parameter d in the range [0, 1/2 [ and for 
positive velocity contrast. 

Note, that the parameter d measures the dissymmetry in the sense that an increasing d implies a 
stronger contrast of the velocity gradients. The parameter d being positive the negative gradient 
of u is smaller in absolute value than the positive gradient. The advantage of this formulation 
compared to the usual profile in the literature as v (x) = sin(27rx) + dsin(47nr) in IITBl is, on one 
hand, to avoid two maxima of the velocity field which does not occur in our pore geometry and, on 
another hand, to control the dissymmetry of the velocity flow by tuning the parameter d without 
change the other parameters as u m and a. The flow uq(x) remains derivable for all \d\ < 0.5 as 
the Stokes flow requires it. Note, however that for d — OA the gradient maximum is about 30. 
This large value implies a large acceleration what is not consistent with the Stokes approximation. 
Therefore, in the bifurcation diagrams of this section only the values of d inferior to 0.4 will be 
relevant for the transport in micropump. 

B. Synchronized transport solutions 

sec:amtw 
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FIG. 19: (a) Continuation of the 1-periodic orbits starting in Fig.[8]at u m = 9 and spanned by the asymmetry 
parameter d. Shown are (red) sq, (black) s m and (green) s a branches, (b) Time evolution of solution 
indicated by circles in panel (a) of the branch s a at d = and d = 0.2. 

fig:perasy 
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FIG. 20: (a) Continuation of the 2-periodic solutions starting in Fig. [8] at u m = 9 and spanned by the 
asymmetry parameter d. Each of the three s a branches splits into two branches when d 7^ leading to 6 
branches. Two of them vanish at a turning point. 

fig:p2cd 
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FIG. 21: Continuation of the synchronized transport solutions starting from Fig. [9] at u m = 9 (7 = 10, 
a = 0.65 and d = 0) and spanned by d. The [black] red branch corresponds to a transport to the [left] right. 
The blue line displays the 2-periodic transport solution to the right. Solid [dashed] lines indicate stable 
[unstable] solutions. Dots designate period doubling bifurcations. 

fig:twu9cd 

One studies the influence of the parameter d on the existence of the synchronized transport 
solutions. Thus, one focusses on the ranges where all periodic solutions are unstable as in Fig. [8} 
The asymmetry of the velocity profile has a poor qualitative influence on the periodic solutions. 



The bifurcation diagram Fig. 19 -a shows that 1-periodic solutions exist for all values of asymmetry 
and there is no bifurcation, for instance they remain all unstable. Note that when d = the 
asymmetric s a branch contains two solutions deduced from the S transformation. Then if d is 
non-zero it splits into two distinct branches. The trajectory of different branches are only slightly 



modified as shows Fig. 19 b. The bifurcation of 2- and 4-periodic solutions reveals a similar 



scenario. The Fig. 20 shows that all 2-periodic solutions are unstable. There is a turning point 
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FIG. 22: Dotted are time evolutions of transport (a) to the left (b) to the right close to the onset of their 
respective existence onset d sn \ and d sn 2. Plain and dashed lines are 2-periodic orbits chosen for the same 



parameters as the transport solution (circles in Fig. 20 1. The different 2-periodic orbits are shifted one 
spatial period and one temporal period. 
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where two branches end but it occurs for d > 0.4 which is not relevant in our framework. Then, 
one may conclude that all periodic solutions remain unstable for relevant values of asymmetry d. 



This conjecture is corroborated by the time-integration of the particle motion Eq. 21 
One expects that the transport solutions found in the symmetric case remains when d is non-zero. 
Indeed the transport solutions already break the reflection symmetry S and thus generically the 
solutions exists, at least, in a neighborhood of d = 0. As for the branch s a , the transport solution 
s t splits into two solutions s t+ and s t _ corresponding respectively to the transport to the right 
(positive x) and to the left (negative x). One starts from the phase point u m = 9, a = 0.65 



and 7 = 10 of the symmetric case (circle in Fig. 10) case and according to the Fig. 10 there is 



one stable transport solution and one unstable. Each transport solution s t+ and s t - is defined by 



respectively Eqs. (35) and (40) of the symmetric case: 



s t+ (t) = s p+ (t)+t (46) 
s t -(t) = s P 4t)-t (47) 

where s p+ and s p _ are 1 -periodic functions. The s t+ and s t - branches are different when d is non- 
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zero. Then, the continuation of the periodic solutions s p+ and s p _ requires two different ODEs: 

Sp+ = 7 (u m (l + aw (s p+ (t) + t, d)) sin(27rt) - 1 - s p+ ) (48) 
s p - = 7 (u m (l + aw (s p -(t) —t,d)) sin(27rt) + 1 - s p _) (49) 

There is four branches starting from both transport solutions of the symmetric case. Stable and 
unstable transport solutions of each direction meet in a saddle-node bifurcation exchanging their 



stability (Fig. [21]) . The left transport ceases at d sn i = 0.302 and the right transport at d sn2 = 0.433 
(Fig. 21). The right transport branch loses its stability in the range [g^l = 0.0133; c^hO. 346] 
via two period doubling bifurcations. Again there is a cascade of bifurcations leading to 
2 n -periodic transport solutions. The transport velocity of these branches remains c = +1. Using 
the time integration, one is found that all these transport to the right become unstable when 
d < du ~ 0.306. Thus d > dn, all dynamics are attracted to a stable transport with c = +1. In 
the narrow range [d sn i, d^\, the transport to the left is no more stable: the simulation shows that 
the dynamics switch from left and right transport as a random walk. Because of the asymmetry, 
generically, the basin of attraction of opposite transport direction are different and then there is 
mean transport direction after a long time integration. The velocity transport is no more locked 
and depends continuously on the parameters, here it decreases with d from 1 to —1. The scenario 
of this transport reversal: transport solution "destroyed" by way of crisis and birth of transport 
solution at a tangent bifurcation is widely described in deterministic inertia ratchet literature for 
instance in 11271 l28l . The continuation diagrams allows to interpret this transport reversal as a 
consequence of the shift of existence domains of right and left transports. This shift is itself due 
to the asymmetry of the pore ratchet (d ^ 0). That the tangent bifurcation occurs in the chaotic 
region of the right transport is a coincidence. There is many other scenarios as we show in the 



next bifurcation diagram Fig. 23 For di < d < d sn i only the transport to the left (c = —1) is 
stable and the simulation shows that it attracts all dynamics. Finally if d < d^, the right transport 
becomes again stable and both transport solutions, to the left and to the right, are stable. Then 
the transport direction depends on the initial conditions. If one considers the initial condition as a 
random variable, generically, one expects a probability different to go to the left or to the right. 
The influence of a lattice of shifted 2-periodic orbits is still relevant in the emergence of the 
transport solution. Among the six 2-periodic branches of continuation diagram in Fig. |20l we 



chose the solution for which two consecutive shifted orbits are locally close. The Fig. 22 shows 



that the transport is in the vicinity of the 2-periodic orbit during one period and switches to 
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FIG. 23: Continuation of the left and right synchronized transport solutions starting in Fig. 21 at d = 0.1 
(a = 0.65 and u m = 9) and spanned by 7. The color, stability and bifurcation codes as in Fig.[2T 
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another shifted 2-periodic orbit as explained for the symmetric case. The transport velocity is 
locked either to +1 or —1. 

Thus by setting the parameter d, i.e. the geometry of the pore, the transport direction may 
be reversed. However in the experiment a switch of transport direction is observed only by 
modifying the pumping frequency, is it possible in our model? The parameter d depending only 
on the pore geometry, it must be fixed. In contrast, a change of the pumping period implies a 
change of the dimensionless drag 7. Thus one expands the transport branches by the parameter 7 



and one starts from the transport solution found in Fig. 21 at d = 0.1. Each branch of transport 
solutions expanded by 7 forms a loop and emerges/ends by a tangent bifurcation (Fig. [23]) as 
for the symmetric case. The left transport exists in the range [yl nL = 7.93;7f ni? = 13.5349] 
and the right transport for 7 G [Y snL = 7.37;7^ ni? = 11.72]. Moreover, for the right transport 
there is two period doubling bifurcations at 7 = 7.96 and 10.64. Both bifurcations are connected 



by a 2-periodic transport solution (Fig. 23). This branch is unstable via two doubling period 
bifurcations (Fig. 23 ). Indeed a cascade of period doubling bifurcations occurs. All these branches 
are transport solutions with the velocity c = +1 as explained in the previous paragraph. For 
example, the time integration at 7 = 10.25 close to the boundary Y R shows a stable transport with 
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FIG. 24: Solutions of synchronized solutions of (upper panel) the branch s t + and (lower panel) the branch 
st_ . Shown is the particle position in the comoving frame [(lower panel) x(t) = —t] (upper panel) x(t) = t 
at [7 = 9.1] 7 = 10.25. The st+ is dynamics is quasi-periodic and s t - is 2-periodic. 
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c = +1. In the comoving frame x(t) = t, the dynamics does not seem to be periodic because a 



cascade of doubling periodic bifurcations occurred for this parameter value (Fig. 24). There is a 



narrow range for which this dynamics is stable. The time integration of the ODE Eq. 21 shows 
that for 7 < Yr — 10.23 and 7 > Yl = 8.01 all periodic branches of 2™ periods are unstable. This 



range is slightly shorter that the instability range of the 2-periodic branch (Fig. 23 ). The left most 
limit 7£ is close to the turning point, then there is an island of stability of the right transport in 
the range [YsnLi 7iJ- A similar scenario occurs for the transport to the left: the left transport with 
c = — 1 is unstable in the range [y l L = 7.94; 7^ = 8.99] due to a period doubling cascade. This 



range is slightly shorter that the range of the 2-periodic solution instability (Fig. 23 ). An example 



of the transport dynamics is showed in Fig. 24 when 7 = 9.1 < 7^. The solution in the comoving 
frame x(t) = — t is a stable 2-periodic transport solution and attracts all the dynamics. Finally, 
close to the leftmost saddle node there is an island of stability of the left transport. 
The bifurcation diagram Fig. [23] allows to answer positively of our question since there is distinct 
ranges of 7 values corresponding to stable transport solutions in different direction. In the 
range [7f ni? ; Y sn R\ on ^y me transport to the right is stable and the time integration shows that the 
dynamics is attracted by this solution. It is an absolute transport to the right. In the same way, in 
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FIG. 25: (a) Time evolution of the particle position at 7 = 8.8 with the parameters of the bifurcation 



diagram in Fig. 23 ( u m = 9, a = 0.65, d = 0.1). (b) Time evolution of the transport velocity c(t) for four 
different initial values. 

fig:simutwg8,8 
the range py£; 7^] it is an absolute transport to the left. 

Between these ranges ([7^; 7f nR ]) there is two stable transports in opposite directions. Depending 
on the initial conditions the dynamics will be attracted by one of these two stable solutions. 
Because of the dissymmetry of the geometry, the attraction basin of each solution is different and 
an effective transport occurs. If one takes the mean velocity over the initial conditions as in ifTBI 
the mean drift velocity varies continuously from 1 to -1 when 7 is decreasing from j^ nR to 7^. 
Such transport reversal was already interpreted as a consequence of bistability in [|2ll28l. 
There is a competition between left and right transport solutions in the range py£; 7^] too. This 
time both transport solutions are unstable. The numerical simulation shows dynamics which 



alternate the transport directions (Fig. 25 -a). Moreover, the switch between the directions lasts 
less than one period, i.e. one does not observe events of quasi-periodic motions. For example, at 
7 = 8.8, the drift is negative and its velocity is in the range c = —0.12 and —0.2 (Fig. |25|b). In 
the remaining range [jj nL ; l\ n j\ only the transport in the left direction exists but it is unstable. 
Then, there is a competition between bound attractors as quasi-periodic solutions and the transport 



solution. Fig. 26 -a displays episodes of transport close to the synchronized transport solution 



following by quasi-periodic motion episodes (tens of periods). One does not observe transport 
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FIG. 26: (a) Time evolution of the particle position at 7 = 7.5 with the parameters of the bifurcation 

9, a = 0.65, d = 0.1). We observe long periods of oscillations in a pore following 



diagram Fig. 23 ( u Tt 

by transport episodes leading to a slow drift to the left, (b) Time evolution of the transport velocity c(t) for 
four different initial values. 
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to the right. Therefore the velocity transport is negative. For example for 7 = 7.5, the drift 



velocity c(t) varies between —0.15 to —0.2 depending on the initial conditions (Fig. 26 -b). The 
mean transport velocity is no more locked and it is presented in literature as a non phase locked 
transport ll2~8ll . However, our analysis shows that synchronized transports exist and generate this 
dynamics. 



To conclude, there is synchronized transport solutions involving 2-periodic orbits as in the 
symmetric case. The existence domain of transport solution is slightly modify for a non-zero 
dissymmetry d. However, the transport may exist in only one direction and may attract the whole 
dynamic ensuring an absolute transport to this direction. As for the experiment it is possible to 
change the transport direction only by varying the frequency of the pumping and not the pore 
lattice geometry. 
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FIG. 27: Particle position at entire times n in the comoving frame (c = +1) for two different values of 7 
close to the tangent bifurcation of synchronized transport. The dynamics displays an irregular staircase for 
which the length increases and diverges when 7 approching the tangent bifurcation. 
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C. Chaotic intermittent transport (7 about 10) 

Beyond the existence domain of synchronized transport there is a complex transport. The 
tangent bifurcations at Y sn R suggests the presence of an intermittent bifurcation in the comoving 
frame. The time integration of the solution in the comoving frame s(t) — t for 7 > Y sn R 
shows long plateaux corresponding to a transport with velocity exactly equal to plus one. These 
episodes of transport are interrupted by a bounded dynamics in the laboratory frame resulting in a 



descending staircase in the comoving frame( Fig. 27 ). The length of plateaux are irregular but are 
statistically longer approaching the tangent bifurcation Y sn R- ^ * s characteristic of an intermittent 



bifurcation and such a staircase is reminiscent of the Devil's staircase. The return map (Fig. 28) 
of the position suggests that is an intermittent bifurcation of type I described by Manneville 
and Pomeau [fTVll : Before the existence of the synchronized transport solution the dynamics is 
nearly tangent and the time spent in the vicinity diverges. When the dynamics escapes from this 
point a chaotic behavior takes place. Such a behavior agrees with the curve shape in the return 



map (Fig. 28) since it is not monotonous and the slope may be large. Thus two close points 
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FIG. 28: (a) Poincare return map of the particle position s n at entire times n of the dynamics close to 



the tangent bifurcation Y sn R °f tne right transport (bifurcation diagram Fig. 23 ). (b) Magnification of the 
dynamics close to the tangent bifurcation (circle in (a) panel) for two values of 7: 7 = 13.535 (red '+' ) and 
7 = 13.53499 (black dots). 
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are after an iteration far away from each other. When one moves away from the critical point 
(7 larger) the episodes of transport and bounded dynamics are equivalent. At 7 about 14.2 a 



stable synchronized transport with c = +1/3, i.e. a 3-periodic transport solution (Fig. 29), at- 
tracts all the dynamics so that it is not possible to observe the end of the complex transport solution. 



The dynamics for 7 smaller than the onset ^\ nL of synchronized left transport solution shows 
switching between transport to the left and bound dynamics as for the intermittent bifurcation 
to the right. However, one does not observe a divergence of duration of the transport solution 
approaching the turning point 7f ni on the contrary the dynamics is similar to the one found in 



Fig. 26 when the (unstable) synchronized solution exists. Moreover the transport episodes are no 
more regular as the synchronized transport. Thus there is not an intermittent bifurcation as the 
turning point as it occurs for the right transport. Note that the time integration does not show the 
existence of a synchronized transport solution and then the continuation method is needed to detect 
the tangent bifurcation which is the birth of the synchronized transport. One found out the onset 
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FIG. 29: Synchronized transport solution with c = 1/3 drift velocity obtained at 7 = 14.2 with the 

9, a = 0.65, d = 0.1). After three temporal periods 



parameters of the bifurcation diagram Fig. [23] ( u m 
the particles moves exactly of one pore length. 
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7 



7.172 of the intermittent transport. The simulation close to this point displays a vanishing 



drift (Fig. 30). The time evolution in the laboratory frame of the particle position evokes again 
a Devil's staircase. The step does not correspond to a periodic state but rather nearly 2-periodic 
solutions. The dynamics becomes chaotic but still bound just before or after the jump between 
plateaux. The duration close to the bound attractor increases when approaching the critical value 
of 7. Therefore we believe the existence of an intermittent bifurcation from a bound attractor 



(Fig. 31 -a). After the transition, there is additionally points in the Poincare map which connect the 
array of attractors allowing the drift (Fig.[3T}b). 

To conclude, complex transport emerges from the chaotic attractors via an intermittent tran- 
sition resulting in a slow drift. Such a dynamics was called in [28 1 aperiodic non phase locked 
transport. As for this latter paper, their existence domain bound includes the synchronized trans- 
port domain. Our analysis allows to distinguish non phase locked transport resulting from the 
competition between synchronized transport solutions as described in previous section and the 
present intermittent transport. Contrary to the literature for which both non phase locked trans- 
ports are not discerned, we are able to determine the birth of synchronized transport. In the next 
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FIG. 30: Particle position at entire times n close to the onset of the intermittent drift (7 ~ 7.1726). Other 

9, a = 0.65, d = 0.1. The drift velocity vanishes 



parameters are as in bifurcation diagram Fig. 23 
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section, we show that such a intermittent bifurcation exists but involving simply the one-periodic 
solutions for large drag. 



D. Quasi-periodic intermittent transport (large 7) 

sec:aint 

For the symmetric case no transport solution may exist when the drag 7 is large (about 100) 
because of the existence of a stable 1-periodic solution as in the bifurcation diagram of Figs. [2} 
However, when the velocity profile is asymmetric, the existence of such a "symmetric" branch 
solution is no more guaranteed. By varying d, one shows, in this section, that one-periodic 
oscillations solutions s , s m and also s a may vanish for large drag 7. 



The diagram bifurcation (Fig. [32]) of the 1-periodic solutions in the axisymmetric case 
(d = 0) is expanded using the velocity contrast a parameter, the drag 7 being setting to 100 and 
the characteristic velocity u m to 9. The 1-periodic solutions so, s m exchange their stability via 



two pitchfork bifurcations connecting the asymmetric branch s a (Fig. 32). It is a similar scenario 
as it occurs in the figures [2] and [7] of the Section IV The s a branch exists in the small ranges 



around a = 0.65 ([a m ; a M ] inset of Fig. 32) and around a = 0.9. 
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FIG. 31: (a) Poincare (x, v) map of the strange attractor before the drift onset (7 = 7.172). The x axis 



is the congruence modulo 1 of the particle position x. Remaining parameters are the same as in Fig. 23 



(b) Poincare (x,v) map of the intermittent drift shown in Fig. 30 after the onset (7 = 7.175). There is 



supplementary points allowing the connection between the shifted attractors. 

fig:poincareg7, 1 7 

Now, let us study the existence of these solutions for an asymmetric velocity profile, i.e., when d 
varies in the range [0;l/2]. It is found that if a is smaller than a critical value a c about 0.507 then 
both branches s an d s m remain for all values of d and their stability is unchanged. For larger 
value a than the critical one a c both branches annihilate in a fold bifurcation at d = df < 0.5. 
According to numerical simulation, the fold bifurcation seems to start at df = 0.5 when a = a c . 
When a is superior to a c , the fold bifurcation corresponds to the vanishing of 1 -periodic solutions 
as shows the bifurcation diagram spanned by d in Fig. [33|-a for a = 0.6 > a c . In this case 



df = 0.1797384. The time integration of the motion equation Eq. 21 for d > df and close to 



the turning point does not show a periodic motion rather a slow drift to the left consisting of 
two distinct evolution phases that take place on distinct timescales: the particle oscillates in the 
vicinity of the vanished periodic solution at the onset and then drift to the next pore to the left. 



The Poincare map in the (s, t) plane (Fig. 34 ) displays a regular descending staircase for different 
values of d: the plateau corresponding to oscillations close to the onset. The plateaux being longer 
when d approaches df. The computation of the drift velocity c as a function of d — df indicates a 



power law dependence: c%{d — df) 1 ! 2 (Fig. 35 -a). This dependence may be interpreted as type-I 
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FIG. 32: Continuation of one-periodic orbits for the symmetric velocity profile vq(x) Eq. [34]by varying a. 
Parameters 7 = 100, u m = 9. Red line is the so branch, balck line is the s m branch and the green line is 
the s a branch. Dashed lines indicate unstable solution. 
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intermittent bifurcation at the turning point described by Pomeau and Manneville ifTTl |22| . To 
corroborate that one is plotted the first return map of the discrete dynamics s n+ i = f(s n ) of the 
particle position at different time t = n. The function / has a parabolic trajectory curve which 
before the bifurcation (d < df) intersects the bisectrix into two fixed-points (periodic orbits s and 
s m ), then for d = d c tangents the bisectrix and finally for d > d c do no intersect the bisectrix. The 
Fig. 36 shows the return map for d > df where the parabolic curve is clearly visible in the intset. 



Such an intermittent bifurcation is already mentioned for deterministic inertia ratchet for current 

reversal [27] with a similar second order evolution equation. However, the intermittent bifurcation 

occurs for smaller friction (smaller 7) and the dynamics far from the dynamics is chaotic as it was 

found in the previous section. This last point contrasts with the regularity of the drift. 

Let us focus on the understanding of this drift regularity. First we study the discrete dynamical 

system 

s n +i = f(s n ) = s n + g(s n ) (50) 
Here, we consider this system as one dimensional problem and then we will discuss the influ- 



ence of the second phase parameter the velocity v n . According to Fig. 36 one has the following 
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FIG. 33: Continuations of one-periodic orbits by varying d for different values of a. Other parameters being 
fixed: 7 = 100, u rn = 9. The color code is as in Fig. [32] 
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properties 

1. g is 1 -periodic, 

2. g is small, 

3. \g'\is not large,. 

4. g is strictly negative 

The three first assumptions are in agreement with the physical assumptions of the drift problem. 
The periodicity is deduced from the spatial periodicity. Due to the large drag the motion is quasi- 
advective and x n+ i — x n should be small then second assumption is checked. Although g is small, 
its derivates could be large. Due to the Stokes flow the acceleration cannot be large and then \g'\ 
is not large. The fourth assumption is related to the fact that there is no 1 -periodic solution. 
If one calls sq the point closest to the bisectrix, if after N iterations sn is the smallest point such 
that s N > s + 1 then one needs either iV or N — 1 iterations to be larger than s + 2 (See 
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FIG. 34: Poincare map in the (s, t) plane of the dynamics at a = 0.6 and different values of d superior to 



the onset df of Fig. 33 The percentage of the relative d variation: (d — df)/dt is indicated below each plot. 
Remaining parameters: 7 = 100, u m = 9. 
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appendix VII B for a proof). We show moreover that this behaviour does not depend on the initial 
condition. That explains the regularity and also the stability since to descend one staircase step, 
N or N — 1 iterations are needed for all dynamics. The simulation corroborates quite well this 
result: for d = 0.18000 the number of iteration varies between 378-379 and for d = 0.17975 
the number of iteration is either 1692 or 1693. Approaching the critical point the variation of 
number of iteration may be great: when d = 0.1797385 the period is approximatively 18000 but 
this number vary of ±500 iterations. That can be explained by the 2-dimensional aspect of the 
dynamics. In fact approaching the critical point the step length vanishes and the presence of the 
velocity in the iteration implies a small error comparing to the ID model. This error leads to an 
important number of iterations since the step length vanishes. 

Therefore, the regularity of the dynamics does not necessarily imply the existence of a periodic 
solution. Even if it exists an unstable periodic orbits the proof of its existence is a complex question 
|6l[8l. Usually, one focuses only on periodic shadowing orbits which are orbit which stay "close" 
to the unstable periodic orbit. In this sense, we show that for large drag the dynamics is intermittent 
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FIG. 35: Logarithmic scale plots of the dependence of d — df on the drift velocity for (a) a = 0.6 and 



(b) a = 0.65. df is the saddle-node onset in the bifurcation diagrams of Fig. 33 Remaining parameters: 
u m = 9, 7 = 100. 
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and quasi-periodic. This regular staircase contrasts with the "Devil's stair case" found in Fig. 27 



and in Fig. 30 



This intermittent bifurcation still exists when the third branch s a is present. Indeed when a m < 



a < olm and d — 0, the branch s a emerges from sq (Fig. 32). In the bifurcation diagram spanned 



by d (Fig. 33 :>) the branch starting from s a splits into two branches which connect the branch s 
and s m in two fold bifurcations. The state s m being close to s a when d = then they connected 
at a turning point df 2 close to zero in the bifurcation diagram spanned by d (Fig. 33?). This 



fold bifurcation is not an intermittent one because of the presence of the second stable branch s a 
(Fig. [33})). In contrast, the fold bifurcation connecting this branch s a and s leads to an intermittent 
drift. When a increases this branch becomes larger and in particular df 2 increases. On the other 
hand, the critical value df% = df decreases. Then, for a = 0.65, we have dfi ~ df 2 — 0.05 
(Fig.|33~}c). Again the intermittent bifurcation occurs for the rightmost turning point, for instance 



df2 in Fig. 33 -c, with the characteristic power law c%(d — df) 1 / 2 (Fig. 35 -b). However the second 



fold bifurcation influences qualitatively the drift: the Poincare map of the dynamics (Fig. 37) 
displays a regular staircase with two different steps. Each plateau corresponds to the oscillations 
close to the fold bifurcation. When a is close to au the periodic states sq and s a collapse and so 
e.g. at a = 0.673 the branch connecting these solutions is very small (Fig.[33|-d). Beyond the value 



48 



Prenrint - March 4. 201 3 




(a) 



(b) 



FIG. 36: Return map of the particle position s n at entire times n of the dynamics close to the intermittent 
bifurcation df for (a) a = 0.6 and (b) a0.65. The '+' symbol indicates the fixed point at df. For the panel 
(b) there is two fixed points according to Fig. [33]panel a = 0.65. The curve passes very close to this point 
but it does not cross the bisectrix. 
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a M the s a branch vanishes and the saddle-node bifurcation occurs at df which increases with a 
(see Figsj33]at a = 0.68 and a = 0.8). 

The numerical simulations shows that the intermittent drift attracts all the dynamics and it 
is stable even far away from the onset. The transport does not change its direction and so the 
drift velocity c remains negative. The function \c(d)\ increases with the asymmetry parameter 
d. However, the maximal value is not larger than 0.1, e.g. for a = 0.65, u m = 9 and d = 0.4, 
one obtains c(0.4) = —0.0654. It is still small comparing to the synchronized transport. It is 
noteworthy that the drift remains for large drag, at least 7 = 500 as shows the loci of bifurcation 
points in (7, df) plane (Fig. 38). Beyond 500 the critical dissymmetry parameter 0^(7) is larger 



than 0.4 and then this result is less relevant for our model. The curves corresponding to a = 0.65 
and a = 0.7 have a sharp corner are due to the competition between both fold bifurcations as it 



occurs in Fig. 33 For smaller drag, the intermittent transport ceases at about 7 = 50 depending 
on the other parameters as the velocity contrast. For a = 0.7 the domain curve displays a turning 
point for 7 ~ 41. This fold bifurcation is due to the fusion/annihilation of two saddle-nodes. The 
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FIG. 37: Poincare map in the (s, t) plane and d = 0.0494 beyond the onset d c i in the bifurcation diagram 



Fig. 33 c. Each plateau is close to a fixed point at the onset. Parameters a = 0.65, u m = 9, 7 = 100. 

fig:gl00um9a065xt 



first saddle-node is coming from the s and s m branches of Fig. 33 The second saddle-node is 



due to 1 -periodic branches starting from d = 0.5. We do not study these branches because they 
occur for non realistic parameter range. Then large asymmetry is not always benefit to the the 
ratchet effect. Anyway, in the experiment, the particles have a large drag at least about 50 then the 
intermittent transport is relevant for this framework. 

One finds transport solution for smaller characteristic velocity too. The transport direction 
remains the same. The existence domain range is a smaller than u m = 9. For instance for u m = 5 



the range of intermittent bifurcation is 7 G [80; 120] (Fig. 39). For 7 = 120 the critical asymmetry 



parameter d is larger than 0.4, thus the lager value are not relevant. 



The intermittent transport is stable and exists in a large range of parameters for u m superior to 
5, large drags about 100 and for large enough velocity contrast about 0.5. The transport is slow 
about 0.1, i.e. at least 10 temporal periods are needed to drift of one pore. A noteworthy aspect of 
this transport solution is that the drift is always to the left direction (negative x). Thus, for a > 
and d > or a < and d > 0, the possible transport intermittent is always to the left direction. 



According to the Fig. 18 the drift direction corresponds to the sense of the slow increase of the 
velocity amplitude. Then, only the parameters of the velocity profile u determines the direction 
of the transport. Since the velocity of the driving flow is related to the pore shape, this latter 
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FIG. 38: Loci (7, df) of intermittent bifurcation at u rn = 9 for three velocity contrasts (fine dashdotted 
line) a = 0.60; (dashed line) 0.65 and (plain line) 0.70. The upper domain bound by each curve is the 
intermittent transport domain. 
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determines the transport direction. This fact contrasts to the synchronized transport solutions for 
which the direction of transport is strongly dependent on frequency, mean velocity, size of particles 
etc... The chaotic intermittent transport found close to the synchronized transport solutions are not 
observed for large 7. Indeed, for large drag it does not exist n-periodic solution with n > 1, then 
this intermittent transport involving complex periodic solutions does not exist. Physically it is due 
to the quasi-advective motion. 

VI. CONCLUSION 

secconc 

In this paper we have studied the process of the deterministic transport of particles in a pore lattice 
filled of viscous liquid. This liquid is pumped back and forth. To explore the particle dynamics 
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FIG. 39: Loci (7, df) of intermittent bifurcation at u m = 5 and a = 0.65. The upper domain bound by the 
curve is the intermittent transport domain. 
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in the parameter space, we have formulated the transport as a bifurcation problem. We focus on 
a generic problem related to such an experiment described in [fT3l . A simple model problem is 
considered in which the only force acting on the particle is the linear Stokes' drag in a quasi-static 
creeping flow. Moreover the problem is one dimensional then the particle cannot rotate and the 
interaction with the pore wall are neglected. Therefore the only nonlinearity is due to the spatial 
flow variations. 

For large particles drag, 7 about 100, their motion is quasi-advective and they get the periodicity 
of the pumping. For a symmetric field flow and a small velocity contrast a there is two types of 
oscillation: one stable s m around the minimum flow velocity and the second one unstable s around 
the maximum. These solutions are invariant by the symmetric transformation S or S m and then 
there are called symmetric solutions. Taking a larger velocity contrast, there is another 1 -periodic 
solution s a emerging from one of symmetric solution branches via a pitchfork bifurcation. This 
new solution breaks spontaneously the symmetry So or S m . In its existence domain, it is always 

52 




Prenrint - March 4. 201 3 



stable. Therefore for all parameters values, one of these three solutions is stable and attracts 
all the dynamics forbidding any transport. It is no more true when one breaks the symmetry 
of the velocity profile, i.e. d > 0, which d measures the difference between the two extrema 
of spatial velocity gradient. All these solutions may end at a fold bifurcation. The numerical 
simulations show that it is an intermittent bifurcation of type I [|22l| leading to a slow and regular 
transport. Close to the onset there is two phases in the transport: many oscillations in the vicinity 
of the onset solution following by a drift to another pore. One shows that the quasi-periodicity 
of the intermittence is due to the large friction (overdamped ratchet) and the creeping flow. Such 
a transport does not seem to be described in the literature of ratchet problems. An interesting 
aspect of this transport is that its existence domain corresponds to a parameter range relevant 
for the micro-pump experiment. The simulation shows that the drift velocity increases with the 
asymmetry velocity profile d. For realistic flow velocity profiles, the drift velocity is smaller than 
0.1, i.e. one pore drift needs 10 temporal periods. The simulation shows that the drift direction is 
determined by the kind of asymmetry: in the sense of small positive velocity gradient and large 
negative gradient. Therefore, for this transport the asymmetry of the pore geometry is the crucial 
parameter which determines its existence and its velocity. 

The experiment displayed sometimes a faster transport and furthermore the direction of the 
drift may change only by tuning either the frequency or the pumping pressure, the pore geometry 
being unchanged. Thus, the intermittent transport cannot explain all these phenomena. By tak- 
ing smaller friction values 7 about 10 — 20, the asymmetric branch s a loses its stability due to a 
cascade of period doubling. All periodic solutions become unstable in a window of parameters 
and the resulting dynamics is quasi-periodic or chaotic and bound. There is an intermittent bi- 
furcation from this strange attractor which leads to slow drift. However this transition does not 
lead to the synchronized transport which have a locked transport velocity to one, i.e. the spatial 
periodicity over the time periodicity. This latter dynamics exists in a range slightly smaller than 
the intermittent transport and emerges/vanishes at a tangent bifurcation. In the comoving frame, 
this bifurcation may lead to an intermittent behavior: regular transport interrupted by chaotic os- 
cillations. However, we show that the synchronized transport solution is due to the presence of 
2-periodic solutions. The spatial and temporal periodicities imply the existence of a lattice of 2- 
periodic shifted orbits. If the extrema of one orbits have a significant difference, two consecutive 
orbits may be locally closed to each other. The trajectory of the synchronized transport is during 
one period close to a 2-periodic orbit and switches to the neighbor orbit. Therefore the drift ve- 
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locity is locked to one. This mechanism possess similarity with heteroclinic connection lfl4l since 
the unstable manifold of a orbit is also the stable manifold of the next 2-periodic solution. This 
latter condition is sufficient for the existence of a synchronized transport. However it differs from 
heteroclinic cycles because the dynamics does not approach asymptotically the periodic orbit. 
There is other transport solution with c = 1/q, q E Z too. For instance one finds drift velocities 
c = 1/2 and c = 1/3. For c = 1/2 we show the transport solution is driven during two period by 
a 3-periodic orbit and then switches to the orbit in the next pore shifted by two temporal periods. 
Therefore the synchronized transport does not emerge from a chaotic dynamics but rather from a 
array of m -periodic orbits. Indeed the chaotic dynamics is a consequence of a cascade of period- 
doubling, -tripling. . . Moreover the synchronized transport solution appears to be quite generic of 
ratchet since only spatial and temporal translation symmetries are required. The other mechanisms 
are due to non-linear phenomena: The existence of the s a branch results from a pitchfork bifurca- 
tion (parity symmetry breaking) and then via a period doubling it is possible to obtain the array of 
2-periodic shifted orbits. 

All synchronized solutions have a range for which they are stable. If it is the unique stable branch 
then it attracts all the dynamics. The branch loses its stability via a doubling period cascade. In 
the comoving frame the dynamics may be quasi-periodic or chaotic but the drift velocity remains 
lock to one. 



As shows Section IV the synchronized transport solution exists from the parity symmetry of 
ratchet and the asymmetric case appears as a perturbation of this case since the transport is related 
to bifurcations from the s a branch. Moreover the solution are quantitatively slightly modified but 
the relevant role of the asymmetry is that the left-right transports have different existence domains. 
Thus the current reversal is a generic phenomenon for slightly asymmetric ratchet because of both 
opposite transports exist in different parameter ranges. There is many scenarios for the current 
reversal depending on relative positions of branches of solutions. 

All these phenomena as synchronized transport or transport reversal are described in deter- 
ministic ratchet problems. However, the transition involving the so-called crisis is evoked for 
different transitions in our case and the emergence of the transport is related only to the stable 
transport. Indeed, in literature only stable solutions are studied even in the detailed analysis of 
ll28ll . In contrast, our study shows that the unstable branches are required to understand the birth 
of transport solution and the transport direction reversal, it allows to point out the relevant role of 
the periodic orbits. 
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Our 1-D model does not take into account the pore boundary in the drag force. However, for 
the experiment design, it may be relevant for instance in the narrow region of the pore. In this 
case the drag force should be smaller than in our model because of the presence of the wall slows 
down particle drag. Then, the synchronized may appear for smaller particles as predicted in our 
model for which a drag force about 10 or 20 corresponds to particle sizes about 5 fxm. Moreover, 
the dependence of the drag coefficient 7 as a function of the position x adds nonlinearity in the 
evolution equation and then in ODE point of view, spontaneous symmetry breaking as period 
doubling may still exist. 

A more realistic simulation requires obviously a 3D simulation for which the particle motion does 
not stay on the axis. Because of particle rotation and particle/pore interactions the problem is 
much more complicated. But this complexity is still consistent with non linear mechanism which 
may lead to the presented transport solutions. 



Even this work is deterministic, it allows to point out the possible role of noise in the emergence 
of synchronized transport. Indeed, if it exists a array of 2-periodic orbits as described in the 



Section IV D, the particle trajectory may switch from a 2-periodic orbit to another close orbit 
thank the noise even if the transport solution does not exist. Because of the asymmetry in the 
phase space of the orbits, there is a preferred direction of the particle. It would be interesting to 
combine non-linear analysis and stochastic analysis for slightly noisy ratchet problems. 



VII. APPENDICES 



A. proof of the boundness of the velocity and acceleration. 



sec:bound 



Let us consider the general form of the ODE Eq. 21 governing the particle motion: 



x + 7X = 7t> (x(£),t) (51) 

where \vq\ is bound by a constant M. 

If x > M then 'x < 0. Therefore it is easy to show that if x(t ) < M at the time t then 
Vt > to : &{t) < M. Moreover if x(to) > M then x decreases. And if x < M and one comes 
back to the previous case. So, for t < t : x(t) < max(x(to), M) 
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FIG. 40: Scheme of the (x n ) and (s n ) maps. The hatched segment is the range of possible values of sj\r-i- 
If Xq — 1 < stv-i as in this scheme then x\ — 1 < sn < %n- Otherwise x% — 1 < sjv < %i — 1- 

One may apply a similar reasonment for the case i < M. So one concludes that \x(t) | is bounded. 
The triangle inequality implies that \x(t)\ < 2 r ymax(\x(to)\, M). 

Let us remark that it is the acceleration may be large large for large drag 7. The reason is that 
there is no assumption of the acceleration of the field flow. Anyway, for a finite 7 the damping 
forbids any velocity particle discontinuity. 



B. Quasi-periodic intermittent transport. 

The intermittent transport for large drag 7 displays a very regular aspect. Let us prove that the 
one-dimensional map 

x n+1 = f(x n ) = x n + g{x n ) (52) 

with the following assumtpions 

1. g is 1 -periodic, 

2. g is small, 

3. \g'\is not large, 

4. g is strictly negative, 

allow to understand this behavior. The 2-dimensional aspect of the dynamics is discussed later. 
We construct the map x n starting from x . Because of g is negative x n is decreasing. We call iV 



the smallest iteration number such that x N < x — 1 (Fig. 40). The assumption 3 implies that / is 
a strictly increasing function. Then x$ — 1 < xn-i implies that xi — 1 < xn- 
Now consider the map s n starting from x\ < sq < xq. We seek the minimal iteration number in 
order such sm < x — 1. Because of / is strictly decreasing then x n+ \ < s n < x n for all n. For 
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instance xn < sn-i < %n-i- The range [xn, xn-i] contains the point x — 1 (Fig. 40) and there 
is two cases depending if sat_i is inferior to x — 1 or not. If s^-i < x — 1 then M = N — 1 
is the wanted iteration number. Otherwise x 1 — 1 < s N < x N < x — I and then M = N (case 



considered for Fig. 40). 



Therefore the map x n modulo 1 passes in the interval ]x\, Xq] every iV — 1 or iV iterations. Because 
xq is unspecified all dynamics have this quasi-periodic behavior. Furthermore, if xq is near the 
tangent bifurcation then ]x±, x ] is small and then the dynamics seems numerically periodic. 

The previous proof does not yield if one considers the two-dimensional discrete dynamical 
system to take into account the velocity. According to the numerical simulations, it acts as a small 
perturbation to the previous process. Approaching the bifurcation point, the intervals [x„,x n+ i] 
vanish then a perturbation of this map even small may have a large influence of the number of 
iterations. It is was it observed when d is very close to d c . However the variation of iterations is 
smaller compared the mean iteration number. 
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